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Abstract. New high-resolution observations reveal that small-scale magnetic flux concentrations 
have a delicate substructure on a spatial scale of 0.1". Their basic structure can be interpreted in 
terms of a magnetic flux sheet or tube that vertically extends through the ambient weak-field or field- 
free atmosphere with which it is in mechanical equilibrium. A more refined interpretation comes 
from new three-dimensional magnetohydrodynamic simulations that are capable of reproducing 
the corrugated shape of magnetic flux concentrations and their signature in the visible continuum. 
Faculae are another manifestation of small-scale magnetic flux concentrations. It is shown that 
the characteristic asymmetric shape of the contrast profile of faculae is an effect of radiative 
transfer across the rarefied atmosphere of the magnetic flux concentration. Also discussed are 
three-dimensional radiation magnetohydrodynamic simulations of the integral layers from the top 
of the convection zone to the mid-chromosphere. They show a highly dynamic chromospheric 
magnetic field, marked by rapidly moving filaments of stronger than average magnetic field that 
form in the compression zone downstream and along propagating shock fronts. The simulations 
confirm the picture of flux concentrations that strongly expand through the photosphere into a more 
homogeneous, space filhng chromospheric field. Future directions in the simulation of small-scale 
magnetic fields are indicated with a few examples from recent reports. 

The second part of these lecture notes is devoted to a few basic properties of magnetic flux tubes 
that can be considered to be an abstraction of the more complicated flux concentrations known from 
observations and numerical simulations. By analytical means we will find that an electrical current 
flows in a sheet at the surface of a flux-tube for which location we also derive the mechanical 
equilibrium condition. The equations for constructing a magnetohydrostatic flux tube embedded in 
a gravitationally stratified atmosphere are derived. It is shown that the expansion of a flux tube with 
height sensibly depends on the difference in the thermal structure between the atmosphere of the flux 
tube and the surrounding atmosphere. Furthermore, we will find that radiative equilibrium produces 
a smaller temperature gradient within the flux tube compared to that in the surrounding atmosphere. 
The condition for interchange stability is derived and it is shown that small-scale magnetic flux 
concentrations are liable to the interchange instability. 
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PART I: THE SIMULATION OF SMALL-SCALE MAGNETIC 

FIELDS 

With "realistic simulations" computational physicists aim at imitating real physical 
processes that occur in nature. In the course of rebuilding nature in the computer, they 
aspire to a deeper understanding of the process under investigation. In some sense 
the opposite approach is taken by computational physicists that aim at separating the 
fundamental physical processes by abstraction from the particulars for obtaining "ideal 
simulations" or an analytical model of the essential physical process. Both strategies are 
needed and are complementary. In the first part of these lecture notes, we mainly focus 
on "realistic simulations" and comparison with observations and provide a more abstract 



background in the second part. 

The term small-scale flux concentration is used here to designate the magnetic field 
that appears in G-band filtergrams as bright tiny objects within and at vortices of 
intergranular lanes. They are also visible in the continuum, where they are called/acM/ar 
points [1], while the structure made up of bright elements is known as the filigree 
In more recent times, the small-scale magnetic field was mostly observed in the G band 
(a technique originally introduced by MuUer flSl) because the molecular band-head of 
CH that constitutes the G band acts as a leverage for the intensity contrast [|4l|5l|6l|71|8l. 
Being located in the blue part of the visible spectrum, this choice also helps improving 
the diffraction limited spatial resolution and the contrast in the continuum. 

Small-scale magnetic flux concentrations are studied for several reasons: 

- Since they make up the small end of a hierarchy of magnetic structures on the solar 
surface, the question arises whether they are "elemental" or whether yet smaller flux 
elements exist. How do they form? Are they a surface phenomenon? What is their 
origin? 

- Near the solar limb they can be identified with faculae, known to critically contribute 
to the solar irradiance variation. 

- They probably play a vital role in the transport of mechanical energy to the outer 
atmosphere, e.g., by guiding and converting magnetoacoustic waves generated by the 
convective motion and granular buffeting or, more directly, by ohmic dissipation. 

STRUCTURE OF SMALL-SCALE MAGNETIC FLUX 
CONCENTRATIONS 

Recent observations of unprecedented spatial resolution with the 1 m Swedish Solar 
Telescope by [jH and ifTOl reveal G-band brightenings in an active region as delicate, 
corrugated ribbons that show structure down to the resolution capability of the instru- 
ment of 0.1", while isolated point-like brightenings exist as well. The structure made 
up of these objects evolves on a shorter than granular time-scale, giving the impression 
of a separate (magnetic) fluid that resists mixing with the granular material. Figure [T] 
shows an example G-band filtergram from the former paper taken in a remnant active 
region plage near disk center. In this region, intergranular lanes are often completely 
filled with magnetic field like in the case marked by the white, horizontal lines in Fig.[T] 
There, and in other similar cases, the magnetic field concentration is framed by a stria- 
tion of bright material, while the central part is dark. Besides examples of ribbon bands. 
Figure [T] shows also an isolated bright point in the lower right corner. 

The graphic to the right hand side of Fig. [T] displays the emergent G-band intensity 
(solid curve) from the cross section marked by the white horizontal lines in the image 
to the left. Also shown are the corresponding magnetographic signal (dashed curve), the 
blue continuum intensity (dotted), and the Ca H-line intensity (dash-dotted). We note 
that the magnetic signal is confined to the gap between the two horizontal white lines. 
The intensities show a two-humped profile. 

This situation reminds of the flux-sheet model and the "bright wall effect" of earlier 
theoretical work. A first quasi-stationary, self-consistent simulation of a small-scale flux 
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FIGURE 1. Left: G-band filtergram showing the ribbon-like shape of magnetic flux concentrations. 
Right: G-band intensity (soHd curve) along the section indicated by the horizontal white lines in the 
image to the left. Also shown are the magnetogram signal, the continuum intensity at 436.4 run, and the 
Ca II H-line intensity down-shifted by 0.2. From Berger et al. |(9l. 

sheet was carried out by Deinzer et al. [[TTlfTll . popularly known as the "KGB-models". 
The basic properties of this model are sketched in Fig. [2} (A first kind of this sketch 
was published by Zwaan [fT3l .) Accordingly, a small-scale flux concentration, either of 
tube or sheet-like shape, is in mechanical equilibrium with the ambient atmosphere, viz., 
the gas plus magnetic pressure of the atmosphere within the tube/sheet balances the gas 
pressure in the ambient (field-free) medium. This calls for a reduced density in the flux 
concentration with respect to the environment, at least in the photospheric part, where the 
radiative heat exchange quickly drives the configuration towards radiative equilibrium, 
hence to a similar temperature at constant geometrical height. The density reduction 
renders the flux tube/sheet atmosphere more transparent, which causes a depression of 
the surface of constant optical depth, as indicated by the surface of Xc = 1 in Fig. [2} In a 
plage or network region, this effect increases the "roughness" of the solar surface, hence 
the effective surface from which radiation can escape, which increases the net radiative 
loss from these areas. 

The graphics to the right hand side of Fig. [2] shows a sketch of the relative intensity 
emerging from this model, viz., the intensity of light propagating in the vertical direction 
as a function of distance from the flux sheet's plane of symmetry. It corresponds to 
the plot on the right hand side of Fig. [T| The similarity between this model and the 
observation is striking. Turning to a narrower flux sheet/tube would result in the merging 
of the two contrast peaks to a single central peak in both, model and observation, i.e., 
to a ribbon band or bright point, respectively. Yet, the striation of the depression wall 
that can be seen in the observation is of course not reproduced by the model, which is 
strictly two-dimensional with translational invariance in lane direction. We will see in the 
section on the physics of faculae that three-dimensional magnetoconvection simulations 
show a rudimentary striation. The physical origin of the striation is still unknown. 

Accordingly, the basic properties of ribbon-like magnetic flux concentrations can be 
understood in terms of a magnetic flux sheet embedded in and in force balance with a 
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FIGURE 2. Sketch of a magnetic flux sheet (left) with corresponding intensity contrast (right), distilled 
from self-consistent numerical MHD simulations. Note that the isothermal surfaces are not exactly parallel 
to the surface of optical depth unity, which gives rise to the particular M-shape of the contrast profile. From 
Steiner L14i . 



more or less field-free ambient medium. This can also be said (replacing the word sheet 
by tube) of the rosette structure visible in other still images of Berger et al. [9J who 
call it "flower-like". Flowers can transmute to pores and vice versa. The striation of their 
bright collar is similar to that seen in ribbon structures. Discarding the striation, the basic 
properties of flowers can well be interpreted in terms of a tube shaped flux concentration 
like the one sketched in Fig.[2]or shown in Fig. 30 



A 2 h sequence of images by Rouppe van der Voort et al. ifTOl with a quality compara- 
ble to Fig.[T]reveals that the shape of the ribbon- like flux concentrations and the striation 
of ribbons and flowers change on a very short time-scale, of the order of the Alfven 
crossing travel time. This suggests that these morphological changes and the striation 
itself are related to the flute instability, which small-scale flux concentrations are liable 
to. For an untwisted axisymmetric flux tube, the radial component of the magnetic field 
at the flux-tube surface must decrease with height, dBr/dz\^ < 0, in order that the flux 
tube is stable against the flute (interchange) instability [15]. While sunspots and pores 
with a magnetic flux in excess of 4> ~ 10^^ Wb do meet this condition, small-scale flux 
concentrations do not fulfill it lfT6l[T7l[T8l . Biinte ||T9l shows that small-scale flux sheets 
too, are flute unstable, and he concludes that filament formation due to the flute insta- 
bihty close to the surface of optical depth unity would ensue. As the flux sheet is bound 
to fall apart because of the flute instability, its debris are again reassembled by the con- 
tinuous advection of strands of magnetic field back to the intergranular lane so that a 
competition between the two effects is expected to take place, which might be at the 
origin of the corrugation of the field concentrations and of the striation of the tube/sheet 
interface with the ambient medium. For a derivation of the flute-instability condition see 
part II of these lecture notes. 

Although the fine structure of small-scale magnetic flux concentrations changes 
on a very short time scale, single flux elements seem to persist over the full du- 
ration of the time sequence of 2 h. They may dissolve or disappear for a short 
period of time, but it seems that the same magnetic flux continually reassembles 
to make them reappear nearby. Latest G-band time sequences obtained with the 
Solar Optical Telescope (SOT) on board of the Japanese space satellite HINODE 



(http ://solarb.msfc.nasa. gov) confirm these findings even for G-band 
bright points of low intensity. This suggest a deep anchoring of at least some of the flux 
elements although numerical simulation seem not to confirm this conjecture. 

As indicated in the sketch of Fig. |2} the magnetic flux concentration is framed by a 
downflow of material, fed by a horizontal flow that impinges on the flux concentration. 
Already the flux-sheet model of Deinzer et al. (TT\ showed a persistent flow of this kind. 
According to these authors it is due to radiative cooling from the depression walls of 
the magnetic flux concentrations (the "hot wall effect") that causes a horizontal pressure 
gradient, which drives the flow. The non- stationary flux-sheet simulations of Steiner 
et al. f2G\ and Leka and Steiner f2V\ showed a similar persistent downflow, which 
becomes faster and narrower , with increasing depth, turning into veritable downflow 
jets beneath the visible surface. While downflows in the periphery of pores have been 
observed earlier [|2T1 l22l l23l and also horizontal motions towards a pore [|24l. only 
very recently such an accelerating downflow has been observationally detected in the 
immediate vicinity of ribbon bands by Langangen et al. [i25il . 

SIMULATION OF SMALL-SCALE MAGNETIC FLUX 
CONCENTRATIONS 

For a realistic simulation of small-scale magnetic flux concentrations in the solar pho- 
tosphere one must solve the system of the magnetohydrodynamic equations including 
continuum, momentum, induction, and energy equation, preferably in three spatial di- 
mensions. Since flux concentrations are predominantly governed by convective motion, 
at least the surface layers of the convection zone must be taken into account in order 
to obtain reliable results for the magnetic field in the photosphere. Hence, the compu- 
tational domain must span the distinct radial section of the sun, where energy transport 
changes from convective to radiative so that radiative transfer must be taken into ac- 
count in the energy equation. An equation of state and opacities appropriate for the solar 
plasma in the region of interest are essential for a realistic simulation. The spatial extent 
of the computational domain is limited by the computational resources and the mini- 
mally required spatial resolution. Since we would like to resolve at least magnetic flux 
concentrations with a size of 100 km we should have a grid constant not surpassing 
about 10 km in the horizontal direction. In the vertical direction the resolution can be 
adjusted according to the variation of the pressure scale height but should also be at least 
!^ 10 km in the photosphere. Considering that a computational grid with 1000^ grid cells 
is definitely the uppermost limit for present day simulations, we obtain a size limit for 
the computational domain of about (10 Mm)-^. 

Fig.|3]shows the typical size of the horizontal extent of a two-dimensional and a three- 
dimensional computational domain as used in the past, in relation to granulation and G- 
band bright points. We should be aware that such a computational box, despite the large 
number of grid cells, encompasses a tiny portion, only, of the solar convection zone. 
This is illustrated in Fig. |4[ where the computational box of the size given to the left, is 
shown in relation to the size of the solar convection zone. Yet, this tiny computational 
box includes the full height range of the photosphere and the mass density varies by 
about a factor of 10^. Even though the box reaches only 1.4 Mm into the convection 



FIGURE 3. G-band filtergram of an area of 32 x 22 arcs of the solar surface. The white line drawn 

across an elongated bright ribbon band corresponds to the width of a two-dimensional simulation domain 
of 2400 km. The white box shows the size of a three-dimensional simulation domain of 4800 x 4800 km. 
Image by K. Mikurda, F. Woger, & O. von der Liihe with the Vacuum Tower Telescope (VTT) at Tenerife. 



zone, it probes a depth range of already about 5 pressure scale heights of the convective 
layers and hydrogen changes from completely neutral to a substantial fraction being 
ionized. 

The ideal MHD-equations can be written in conservative form as: 



dt 



(1) 



where the vector of conserved variables, U, the source vector, S, due to gravity and 
radiation, and the flux tensor, are 



U= (p,pv,B,£) , S = (0,pg,0,pg-v + ^rad) , 
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and where the dyadic tensor product of two vectors a and b is the tensor ab = C with 
elements Cmn — ambn- v and B are vectors of the velocity and magnetic fields. jU is the 
magnetic permeability, which can be taken to be jU = jUq = 4;r • 10~^ [V s m~^] for 
the solar plasma. I is the identity matrix and a • b = Y^k^k^k the scalar product of the two 
vectors a and b. 

The total energy E is given by 
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FIGURE 4. Typical size of a three-dimensional computational box (left) on scale with the convection 
zone boundaries of the sun (right). 



where e is the thermal energy per unit mass. The additional solenoidality constraint, 

V-B = 0, (5) 

must also be fulfilled. The MHD equations must be closed by an equation of state which 
gives the gas pressure, p, as a function of the density and the thermal energy per unit 
mass 

p = p{p,£), (6) 

usually available to the program in tabulated form. 
The radiative source term is given by 

qmd = 4;rp jKy{jy-By)dV, (7) 

where Jy is the mean intensity and By the Planck function. Assuming strict local 
thermodynamic equilibrium (LTE), the intensity is given by the formal solution of the 
radiative transfer equation: 

Mr) = ^flv{r,n)da, /v(r,n) = /voe"^^° + ^^'°5v(Tv)e-^MTv , (8) 

where Tyo is the optical depth at frequency v from the boundary to the location r along 
the direction n. dr = Kds, where K is the total opacity and s the spatial distance along a 
line through location r in direction n. Formally, this can be written more compact as 

7v(r)=Av(r,r')5v(rO+Gv(r), (9) 

where Ay is the integral operator which adds the intensities at r caused by emission at all 
the points r' in the considered computational domain, and where Gy is the transmitted 



mean intensity due to the incident radiation field into this domain. In practice the 
frequency integration is either replaced by using frequency mean quantities (Rosseland 
mean opacity) or it is approximated by a method of multiple frequency bands [|26l . For 
a detailed description of the latter method see ll27ll and [|28ll . 

In practice it is not the ideal MHD-equations that are solved but rather some kind of a 
viscous and resistive form of the equations with flux tensor 
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whereR = vp[(Vv) + (Vv)^-(2/3)(V-v)I] is the viscous stress tensor, t] = (v/Prm) = 
l/(/ia) the magnetic diffusivity with o being the electric conductivity, and T] (j x B) = 
(77//i)(VxB)xB.R-vis the tensor product VmkVk- Prm is the magnetic Prandtl number, 
qturb ^ turbulent diffusive heat flux, which would typically be proportional to the 
entropy gradient: q'"*^^ = — (l/Pr)vprV5, where Pr is the Prandtl number. 

Typically, v is not taken to be the molecular viscosity coefficient but rather some 
turbulent value that takes care of the dissipative processes that cannot be resolved by the 
computational grid. Such subgrid-scale viscosities should only act where strong velocity 
gradients lead to turbulence. Therefore, they typically depend on velocity gradients like 
in the Smagorinsky-type of turbulent viscosity where 
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and where c is a free parameter. This parameter is normally chosen as small as possible 
just in order to keep the numerical integration stable and smooth, but otherwise having 
no effect on large scales. 

Advanced numerical (high resolution) schemes feature an inherent "dissipation" (lim- 
iters) that acts similar to the explicit dissipative terms shown in the above flux tensor, 
Eq. ( [To] ). This artificial viscosity is made as small as possible but just large enough in or- 
der to keep the numerical scheme stable and keep gas pressure and total energy positive. 
One then only has to program the ideal equations. Of course, in this case it is difficult to 
quote the actual Reynolds and Prandtl number of the simulation because these numbers 
change then form grid cell to grid cell depending on the flow. Therefore, for some appli- 
cations, it might be preferable to explicitly include the dissipative terms in the equations 
using constant dissipation coefficients, which then allows for well defined dimensionless 
numbers. However one integrates the ideal equations on a discrete computational grid. 
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FIGURE 5. Typical boundary conditions for thermal variables and velocities. 



one is locked with a discretization error that normally assumes a form similar to the 
dissipative terms in the non-ideal equations. More details on computational methods for 
astrophysical fluid flow can be found in LeVeque [|29ll . Other very useful text books on 
the numerical solution of the gasdynamics equations are [|30ll3n[32ll . For corresponding 
solution methods of the MHD equations see, e.g., |l33l[34l[35l[36l|371|38l and references 
therein. 

Fig. |5] displays typical boundary conditions for thermal variables and velocity. Note- 
worthy are the conditions at the bottom boundary, where the plasma should freely flow 
in and out of the computational domain (subject to the constraint of mass conservation) 
since this boundary is located within the convection zone. Inflowing material has a given 
specific entropy that determines the effective temperature of the radiation leaving the do- 
main at the top, while the outflowing material carries the entropy it instantly has. This 
boundary condition reflects that inflowing material originates from the deep convection 
zone with a close to adiabatic stratification of a fixed entropy. 
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FIGURE 6. Two different realizations for the boundary conditions of the magnetic field. The condition 
of the right hand panel allows for the transport of horizontal magnetic field of a given strength into the 
box. 

Fig. |6] displays two different realizations (both presently in use) for the boundary 
conditions of the magnetic field. While the magnetic field in the condition of the left 
panel is forced to become vertical at the top and bottom boundary, the condition of 
the right hand panel allows for the transport of horizontal magnetic field of a given 
strength into the box. More details about this condition are given in the chapter on future 
directions. 



Results from multidimensional simulations 



New results from realistic three-dimensional simulations on the formation, dynamics 
and structure of small-scale magnetic flux concentrations have recently been published 
in a series of papers by Schiissler and collaborators. Vogler et al. [|39l simulate magne- 
toconvection in a box encompassing an area on the solar surface of 6 x 6 Mm^ with a 
height extension of 1400 km, reaching from the temperature minimum to 800 km below 
the surface of optical depth unity. Although this is only 0.4% of the convection zone 
depth, the box still includes the entire transition from almost completely convective to 
mainly radiative energy transfer and the transition from the regime where the flux con- 
centration is dominated by the convective plasma flow to layers where the magnetic 
energy density of the flux concentrations by far surpasses the thermal energy density. 
The bottom boundary in this and similar simulations is open in the sense that plasma 
can freely flow in and out of the computational domain, subject to the condition of mass 
conservation. Inflowing material has a given specific entropy that determines the effec- 
tive temperature of the radiation leaving the domain at the top, while the outflowing 
material carries the entropy it instantly has. 

Figure |7] shows a snapshot from this simulation: To the left the emergent mean 
intensity, to the right the vertical magnetic field strength at a constant height, viz., at the 
horizontally averaged geometrical height of optical depth unity. (I would like to caution 
that this magnetic map is not what would be seen with a magnetograph, irrespective 
of its spatial resolution, because it refers to a plane parallel section, which is not what 
is sampled with a magnetographic or polarimetric mapping.) The strong magnetic field 
in intergranular lanes is manifest in a corresponding signal in the emergent intensity 
very much like the observations shown in the previous sections (Figs. [T] and [3]). Also, the 
intensity signal shows the same corrugated and knotted ribbon structure that is observed, 
and sometimes there appear also broader ribbon structures with a dark central core, 
like the one marked in Fig. [T] The characteristic striation of the latter case however, 
is absent in the simulation, possibly because the flute instability is suppressed on very 
small spatial scales due to lack of sufficient resolution of the simulation. In the central 
part of the snapshot, a micro pore or magnetic knot has formed. 

A comparison of the average gas plus magnetic pressure as a function of height at 
locations of magnetic flux concentrations with the run of the average gas pressure in 
weak-field regions reveals that the two are almost identical, proving that even in this 
dynamic regime, the thin flux tube approximation is very well satisfied [39J. This result 
confirms that the model discussed in the previous section and sketched in Fig.|2]is indeed 
an acceptable first approximation to the real situation. 

Simulations are not just carried out for the sake of reproducing observed quantities. 
Once good agreement with all kind of observations exists, simulations allow with some 
confidence to inform about regions not directly accessible to observations, for exam- 
ple about the magnetic structure in subsurface layers. In this respect the simulations 
of Vogler et al. [39] show that often flux concentrations that have formed at the sur- 
face disperse again in shallow depths. This behaviour was also found by S chaff enberger 
et al. pTI in their simulation with an entirely different code and further by Stein and 
Nordlund [|42l . A vertical section through a three-dimensional simulation domain of 
Schaff enberger et al. bMjl, where two such shallow flux concentrations have formed, is 
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FIGURE 7. Simulation snapshot. Left: Frequency integrated intensity. Right: Vertical magnetic flux 
density at constant average geometrical height of optical depth unity. The mean flux density is 0.02 T. 
From Vogler et al. 



shown in Fig. 14 The superficial nature of magnetic flux concentrations in the simu- 
lations, however, is difficult to reconcile with the observation that many flux elements 
seem to persist over a long time period. 

Although the turbulent nature of convection requires simulation in three-dimensional 
space, it can be instructive to carry out less costly two-dimensional calculations that 
allow for high spatial resolution. For example, a two-dimensional approach can be 
justified for the simulation of elongated sheets of magnetic flux concentrations as they 
can be found in intergranular lanes of active regions, e.g., the G-band bright string that 
extends perpendicular to the 2-D domain indicated in Fig. [3] 

Fig. [8] (top) shows a time instant of a two-dimensional simulation of a magnetic flux 
sheet in interaction with convective motion EOl . A granule has grown to the right of 
the flux sheet and pushes it to the left. Such "granular buffeting" can produce waves 
propagating along, within, and in the environment of a magnetic flux sheet. Besides a 
swaying motion of the flux sheet that excites transverse tube waves, the simulation also 
revealed that often shock waves propagate in the vertical direction in the tenuous upper 
photospheric part within the flux sheet . Such a shock wave is visible in Fig. [8] (top) in a 
height range between 300 to 400 km slightly left of the central part of the flux sheet. 

From simulations as the ones shown in Figs. |7] or [8] (top), synthetic observable 
quantities like intensity contrasts, or spectral lines, or polarimetric quantities can be 
computed and directly compared with real observations. For example. Fig. [8] (bottom) 
displays Stokes V profiles computed for vertical lines of sight distributed over the 
horizontal interval between 600 and 1380 km of the simulation snapshot. The shock 
wave produces the pathological Stokes V profiles belonging to .x; = 900 and 920 km. 
There is a systematic change in the area asymmetry of the Stokes V profiles across the 
flux sheet. While the central part of the flux sheet produces profiles with a negative 
area asymmetry, the peripheral profiles are positive. This behaviour can be understood. 
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FIGURE 8. Top: Snapshot of a two-dimensional simulation of a magnetic flux sheet in interaction with 
convective motion. Color coding indicates temperature in K, white arrows the velocity. Bottom: Stokes V 
profiles that emanate from vertical lines of sight distributed over the horizontal interval between 600 and 
1380 km of the simulation snapshot. Positive area asymmetries result in the two regions labeled with "+", 
while the lines of sight in the middle region labeled with "— " contribute negative area asymmetries. The 
actual profiles emanating from each line of sight is also shown, where the mean amplitude, the amplitude 
asymmetry, the area asymmetry (all in %), and the zero-crossing shift in mA are printed (from top to 
bottom) on the right hand side of each panel. In addition, positive and negative area asymmetries are 
annotated in the right margin with a plus and a minus sign, respectively. From Steiner BOl . 



noting that the flux sheet is most of the time framed by a strong downflow of material 

As mentioned in the first chapter, such a downflow is indeed visible in high resolu- 
tion Dopplergrams of small-scale magnetic flux concentrations [|25| . Since polarimetric 
measurements have not yet reached a similarly good quality in spatial resolution, the 
corresponding polarimetric signal of positive Stokes V area asymmetry in the periphery 
of flux concentrations has not yet been observationally confirmed, but I am confident 
that it soon will be. 

Other results from the computation of synthetic Stokes profiles from three- 
dimensional magnetohydrodynamic simulations have recently been published by 
[143J and [|44|. Further recent three-dimensional magnetohydrodynamic simulations of 
small-scale solar magnetic fields include the works [42J and [|45]| . 

THE PHYSICS OF FACULAE 

With growing distance from disk center, small-scale magnetic flux concentrations grow 
in contrast against the quiet Sun background and become apparent as solar faculae close 
to the limb. Ensembles of faculae form plage and network faculae that are as conspicuous 
features of the white light solar disk, as are sunspots. It is therefore not surprising that 
they play a key role in the solar irradiance variation over a solar cycle and on shorter 
time scales [46i,l47ll48l . Measurements of the center to limb variation of the continuum 
contrast of faculae are diverse, however, as the contrast is not only a function of the 
heliocentric angle, ji = cos 6, but also of facular size, magnetic field strength, spatial 
resolution, etc., and as measurements are prone to selection effects. While many earlier 
measurements report a contrast maximum around jU ~ 0.2 . . . 0.4 with a decline towards 
the limb, latest measurements [|49l [50l [5T1 l52l l53l l54l point rather to a monotonicaUy 
increasing or at most mildly decreasing contrast out to the limb. 

The standard facula model [55 J, again consists of a magnetic flux concentration em- 
bedded in and in mechanical equilibrium with a weak-field or field-free environment as 
is sketched in Fig, [2} When approaching the limb, the limb side of the bright depression 
wall becomes ever more visible and ever more perpendicular to the line of sight, which 
increases its brightness compared to the limb darkened environment. At the extreme 
limb, obscuration by the centerward rim of the depression starts to take place, which 
decreases the size and possibly the contrast of the visible limb-side wall. 

Recently, Lites et al. Il56l and Hirzberger and Wiehr [f57l have obtained excellent 
images of faculae with the 1 m Swedish Solar Telescope. Figure |9] (taken from [|57l ) 
shows on the left hand side network faculae at a heliocentric angle of /i = 0.48 in 
the continuum at 587.5 nm. The solar limb is located towards the right hand side. It 
is clearly visible from this image that faculae are in reality partially brightened granules 
with an exceptionally dark and wide intergranular lane ("dark facular lane") on the disk- 
center side of the contrast enhancement, which is also the location of the magnetic flux 
concentration. The right half of Fig. [9] shows the string of faculae within the white box 
of the image to the left, aligned according to the position of the dark lane. Also shown is 
the mean contrast profile, averaged over the alignment. Similar contrast profiles of single 
faculae are shown by Lites et al. [.56,1 . Such contrast profiles pose now a new constraint 



FIGURE 9. Left: Network faculae at a heliocentric angle of = 0.48 in the continuum at 587.5 nm. The 
solar limb is to the right. Right: Faculae within the white box of the image to the left aligned according to 
the position of the dark lane, together with the mean spatial scan through the aligned faculae. The dotted 
curve refers to profiles for which the maximum intensity exceeds 1.3. From Hirzberger and Wiehr l(57l . 



that any model of faculae must satisfy. 

Magnetoconvective simulations as the one discussed in the previous section indeed 
show facular-like contrast enhancements when computing the emergent intensity along 
lines of sight that are inclined to the vertical direction for mimicking limb observation. 
Such tilted three-dimensional simulation boxes are shown in the papers by Keller et al. 
[I58l . Carlsson et al. [|59||, and De Pontieu et al. [Wj. Keller et al. [58J also show 
the contrast profile of two isolated "faculae", which however have a more symmetric 
shape, rather than the observed characteristic steep increase on the disk-center side 
with the gentle decrease towards the limb. Also they obtain a maximum contrast of 
2, far exceeding the observed value of about 1.3. It is not clear what the reason for 
this discrepancy might be. Interestingly, the old "KGB-model" [fT2l [6Tll as well as 
the two-dimensional, non- stationary simulation of Steiner [l62ll do nicely reproduce the 
asymmetric shape and the dark lane. 

Another conspicuous property of faculae that high-resolution images reveal is that 
they are not uniformly bright but show a striation not unlike to and possibly in con- 
nection with the one seen in G-band ribbons at disk center. While this feature can- 
not be reproduced in a two-dimensional model, it must be part of a satisfactory three- 
dimensional model. But so far 3-D simulations show only a rudimentary striation. This 
finding, rather surprisingly, indicates that the effective spatial resolution of present-day 
three-dimensional simulations is inferior to the spatial resolution of best current obser- 
vations. 

In an attempt to better understand the basic properties of faculae, Steiner [|62ll con- 
siders the ideal model of a magnetohydrostatic flux sheet embedded in a plane parallel 
standard solar atmosphere. For the construction of this model the flux-sheet atmosphere 
is first taken to be identical to the atmosphere of the ambient medium but shifted in the 
downward direction by the amount of the "Wilson depression" (the depression of the 
surface of continuum optical depth unity at the location of the flux concentration). The 
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FIGURE 10. Left: a) Magnetic flux concentration (blue, verticafly oriented lines of force, see internet 
version for colours) with surfaces of optical depth T = 1 for vertical lines of sight (thick/red curve) and 
lines of sight inclined by = 60° to the vertical direction (thin/red curve), b) Corresponding contrast 
curves. All values of the light red curve left of the black dot originate from lines of sight left of the one 
indicated in panel a). Right: a) Surfaces of optical depth T = 1 and 5 (thick/red) for lines of sight inclined 
by 50° to the vertical, together with isotherms, b) Contrast profile. The region of negative contrast is 
bounded by the two oblique lines of sight indicated in panel a). Adapted from Steiner li62J . 



shifting results in a flux-tube atmosphere that is less dense and cooler than the ambient 
medium at a fixed geometrical height. In the photospheric part of the flux concentration, 
thermal equilibrium with the ambient medium is then enforced. Denoting with index i 
the flux-sheet atmosphere and with e the ambient atmosphere and with W the depth of 
the "Wilson depression", we therefore have 

{Te{z + W) for T,>1 

where Xc is the optical depth in the visible continuum and p,(z) < pe{z) Vz. The lower 
left panel of Fig. [10] shows this configuration together with two surfaces of optical 
depth unity, one for vertical lines of sight (disk center), the other for lines of sight 
running from the top right to the bottom left under an angle of = 60° to the vertical, 
like the one indicated in the figure. The upper left panel shows the corresponding 
continuum enhancement for disk center (double humped profile) and 9 = 60°. Of the 
curve belonging to 6 = 60°, all values left of the black dot belong to lines of sight left of 
the one indicated in the lower panel. This means that the contrast enhancement extends 
far beyond the depression proper in the limbward direction, exactly as is observed. The 



reason for this behaviour is explained with the help of Fig. 1 1 as follows. 

A material parcel located in the solar atmosphere and lateral to the flux sheet "sees" 
a more transparent atmosphere in the direction toward the flux sheet compared to a 
direction under equal zenith angle but pointing away from it because of the rarefied 
flux-sheet atmosphere. Correspondingly, from a wide area surrounding the magnetic 
flux sheet or flux tube, radiation escapes more easily in the direction towards the flux 
sheet so that a single flux sheet/tube impacts the radiative escape in a cross-sectional 
area that is much wider than the magnetic field concentration proper. 



The right hand side of Fig. 10 shows a similar situation as to the left but for a flux sheet 
that is twice as wide. The continuum contrast for lines of sight inclined by = 50° to 
the vertical is shown in the top panel. It can be seen that a dark lane of negative contrast 
occurs on the disk-center side of the facula. It arises from the low temperature gradient 
of the flux-sheet atmosphere in the height range of Tc = 1 and its downshift relative to 
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FIGURE 11. Photons preferentially escape along the line of sight to the right that traverses the magnetic 
flux sheet/tube in comparison to the line of sight to the left under equal zenith angle, because of the rarified 
(less opaque) atmosphere in the flux sheet/tube. From Steiner I14L 

the external atmosphere in combination with the inclined lines of sight. One could say 
that the dark lane in this case is an expression of the cool "bottom" of the magnetic flux 
sheet. 

It is remarkable that this basic, energetically not self-consistent model is capable of 
producing both, the facular dark lane and the asymmetrically shaped contrast curve of 
the facula, with realistic contrast values. The results of this basic hydrostatic model carry 
over to a fully self-consistent model of a magnetic flux sheet in dynamic interaction with 
non- stationary convective motion [|62l . In this case the facular lane becomes broader and 
darker. 

It follows from these insights that a facula is not to be identified with bright plasma 
that sticks, as the name may imply, like a torch out of the solar surface and as the 
"hillock model" of Schatten et al. [|63l suggests. Rather is it the manifestation of pho- 
tospheric granulation, seen across a magnetic flux concentration — granulation that ap- 
pears brighter than normal in the form of so called "facular granules". Interestingly, 
already Chevalier [64] wondered: "La granulation que Ton voit autour des taches plus 
eclatante que sur les autres parties est-elle la granulation des facules ou celle de la pho- 
tosphere vue a travers les facules ?" and Ten Bruggencate [65] noted that "Sie [Pho- 
to spharengranulen und Fackelgranulen] unterscheiden sich nicht durch ihre mittlere 
Grosse, wohl aber durch den Kontrast gegeniiber der Umgebung." 

If this is true, one expects facular granules to show the same dynamic phenomena 
like regular granulation. Indeed, this is confirmed in a comparison of observations with 
three-dimensional simulations by De Pontieu et al. [60]. Fig. [12] shows on the left hand 
side an observed time sequence of a facula in which a dark band gradually moves from 
the limb side toward the disk center (down) and seemingly sweeps over and "erases" the 
facula temporarily. The panel on the right hand side shows the same phenomenon in a 
time sequence from a three-dimensional simulation. At t = s a dark band is just above 
the center of the figure. It moves toward the disk center (downward in the figure), and 
the faculae are reduced to a weak band at t = 90 s. 

Examination of the simulation sequence reveals that dark bands are a consequence 
of the evolution of granules. Often granules show a dark lane that usually introduces 
fragmentation of the granule. The smaller fragment often dissolves (collapses) in which 
case the dark lane disappears with the collapsing small fragment in the intergranular 
space. Exactly this process can lead to the dark band phenomenon, when a granular 
dark lane is swept towards the facular magnetic flux concentration. Since the facular 
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FIGURE 12. Le/f, observation: A dark band gradually moves toward the disk center (down) and 
seemingly sweeps over and "erases" the facula temporarily. Right, simulation: At t = s a dark band is 
just above the center of the figure. It moves toward the disk center (downward), and the facula is reduced 
to a weak band at t = 90 s. From De Pontieu et al. fEU]. 



brightening is seen in the disk-center facing side of granules, only granular lanes that 
are advected in the direction of disk center lead to facular dimming. This example nicely 
demonstrates how regular granular dynamics when seen across the facular magnetic field 
can lead to genuine facular phenomena. 

Despite the major progress that we have achieved in understanding the physics of 
faculae over the past few years, open questions remain. These concern 

- a comprehensive model of the center to limb variation of the brightness of faculae 
including dependence on size, magnetic flux, flux density, color, etc., 

- a quantitative agreement between simulation and observation with respect to mea- 
surements in the infrared and with respect to the observed geometrical displacement 
between line core and continuum filtergrams of faculae, 

- the physical origin of the striation, 

- a quantitative evaluation of the heat leakage caused by faculae, or 

- the role of faculae in guiding magnetoacoustic waves into the chromosphere. 




3-D MHD SIMULATION FROM THE CONVECTION ZONE TO 

THE CHROMOSPHERE 

For investigating the connection between photospheric small scale magnetic fields and 
the chromosphere, S chaff enberger et al. [41 J have extended the three-dimensional radia- 
tion hydrodynamics code CO^BOLEj^to magnetohydrodynamics for studying magneto- 
convective processes in a three-dimensional environment that encompasses the integral 
layers from the top of the convection zone to the mid chromosphere. The code is based 
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FIGURE 13. Horizontal sections through the three-dimensional computational domain. Color coding 
displays log|B| with individual scaling for each panel (see internet version for colours). Left: Bottom 
layer at a depth of 1210 km. Middle: Layer 60 km above optical depth = 1. Right: Top, chromospheric 
layer in a height of 1300 km. White arrows indicate the horizontal velocity on a common scaling. Longest 
arrows in the panels from left to right correspond to 4.5, 8.8, and 25.2 km/s, respectively. Rightmost: 
Emergent visible continuum intensity. From Schaffenberger et al. ||4T|. 



on a finite volume scheme, where fluxes are computed with an approximate Riemann- 
solver [|29l l30ll for automatic shock capturing. For the advection of the magnetic field 
components, a constrained transport scheme is used. 

The three-dimensional computational domain extends from 1400 km below the sur- 
face of optical depth unity to 1400 km above it and it has a horizontal dimension of 
4800 X 4800 km. The simulation starts with a homogeneous, vertical, unipolar magnetic 
field of a flux density of 1 mT superposed on a previously computed, relaxed model 
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FIGURE 14. Snapshot of a vertical section through the three-dimensional computational domain, 
showing log|B| (color coded) and velocity vectors projected on the vertical plane (white arrows). The 
b/w dashed curve shows optical depth unity and the dot-dashed and solid black contours j3 = 1 and 100, 
respectively. See internet version for colours. From Schaffenberger et al. L41J . 



of thermal convection. This low flux density is representative for magnetoconvection in 
a very quiet network-cell interior. The magnetic field is constrained to have vanishing 
horizontal components at the top and bottom boundary, but lines of force can freely 
move in the horizontal direction, allowing for flux concentrations to extend right to the 
boundaries. Because of the top boundary being located at mid-chromospheric heights, 
the magnetic field is allowed to freely expand with height through the photospheric lay- 
ers into the more or less homogeneous chromospheric field. 

Figure [13] shows the logarithmic absolute magnetic flux density in three horizontal 
sections through the computational domain at a given time instant, together with the 
emergent Rosseland mean intensity. The magnetic field in the chromospheric part is 
marked by strong dynamics with a continuous rearrangement of magnetic flux on a time 
scale of less than 1 min, much shorter than in the photosphere or in the convection- 
zone layers. There, the field has a strength between 0.2 and 4 mT in the snapshot of 
Fig. [13} which is typical for the whole time series. Different from the surface magnetic 
field, it is more homogeneous and practically fills the entire space so that the magnetic 
filling factor in the top layer is close to unity. There seems to be no spatial correlation 
between chromospheric flux concentrations and the small-scale field concentrations in 
the photosphere. 

Comparing the flux density of the panel corresponding to z = 60 km with the emergent 
intensity, one readily sees that the magnetic field is concentrated in intergranular lanes 
and at lane vertices. However, the field concentrations do not manifest a corresponding 
intensity signal like in Fig. [7] (left). This is because the magnetic flux is too weak to 
form a significant "Wilson depression" (as can be seen from Fig.[T4|) so that no radiative 
channeling effect (see radiative equilibrium solution in part II) takes place. 

Figure [14] shows the logarithm of the absolute field strength through a vertical section 
of the computational domain. Overplotted are white arrows indicating the velocity field. 
The b/w dashed curve corresponds to the optical depth unity for vertical lines of sight. 
Contours of the ratio of thermal to magnetic pressure, j8, for j8 = 1 (dot-dashed) and 
j8 = 100 (solid) are also shown. Magnetoacoustic waves that form transient filaments of 
stronger than average magnetic field are a ubiquitous phenomenon in the chromosphere 



and are also present in the snapshot of Fig. 14 e.g., along the contour of /3 = 1 near 
X — 1200 km and x = 2500 km. They form in the compression zone downstream and 
along propagating shock fronts. These magnetic filaments that have a field strength 
rarely exceeding 4 mT, rapidly move with the shock fronts and quickly form and dissolve 
with them. The surface of /3 = 1 separates the region of highly dynamic magnetic fields 
around and above it from the more slowly evolving field of high beta plasma below it. 
This surface is located at approximately 1000 km but it is corrugated and its local height 
strongly varies in time. 

A very common phenomenon in this simulation is the formation of a 'magnetic 
canopy field' that extends in a more or less horizontal direction over expanding granules 
and between photospheric flux concentrations. The formation of such canopy fields 
proceeds by the action of the expanding flow above granule centres. This flow transports 
'shells' of horizontal magnetic field to the upper photosphere and lower chromosphere, 
where shells of different field directions may be pushed close together, leading to a 
complicated network of current sheets in a height range from approximately 400 to 
900 km. 
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FIGURE 15. Logarithmic current density, log|y|, in a vertical cross section (top panel) and in four 
horizontal cross sections (bottom panels) in a depth of 1 1 80 km below, and at heights of 90 km, 610 km, 
and 1310 km above the mean surface of optical depth unity from left to right, respectively. The arrows in 
the top panel indicate the magnetic field strength and direction. The dashed line indicates the position of 
the vertical section, [j] = 3 x 10^ A/m^. From Wedemeyer-Bohm et al. Il66l . 



This network can be seen in Fig. 15 (top), which shows, for a typical snapshot of the 
simulation, the logarithmic current density, log|j|, together with arrows indicating the 
magnetic field strength and direction. Figure 15 (bottom) shows from left to right log \j\ 
in four horizontal cross sections in a depth of 1180 km below, and at heights of 90 km, 
610 km, and 1310 km above the mean surface of optical depth unity. Higher up in the 
chromosphere (rightmost panel), thin current sheets form along shock fronts, e.g., in the 
lower left corner near x=lA Mm. 

Using molecular values for the electrical conductivity, Wedemeyer-Bohm et al. [[66ll 
derive an energy flux of 5 to 50 W m^^ into the chromosphere caused by ohmic dissipa- 
tion of these current sheets. This value is about two orders of magnitude short of being 
relevant for chromospheric heating. On the other hand, the employed molecular values 
for the conductivity might be orders of magnitude too high for to be compatible with 
the effective electrical conductivity of the numerical scheme determined by the inher- 



ent artificial diffusion. Therefore, the ohmic heat flux might be conceivably two orders 
of magnitude larger than suggested by this rough estimate, so that magnetic heating by 
ohmic dissipation must be seriously taken into account. More advanced simulations, 
taking explicit ohmic diffusion into account will clarify this issue. The origin of current 
sheets and their microscopic nature is treated in the first chapter of the second part of 
these lecture notes. 

FUTURE DIRECTIONS 

Continuously increasing power of computational facilities together with steadily im- 
proving computational methods, expand the opportunity for numerical simulations. On 
the one hand, more detailed physics can be included in the simulation, on the other hand 
either the computational domain or the spatial and temporal resolution can be increased. 
In most simulations, especially when the computational domain encompasses only a 
small piece of a star, boundary conditions play an important role. They convey informa- 
tion on the outside world to the physical domain of the simulation. But this outside world 
is often poorly known. In order to acquire experience and intuition with respect to the 
influence of different types of boundary condition on the solution, one can implement 
and run various realizations of boundary conditions, which however also requires addi- 
tional resources in computer power and time. Also the initial condition may critically 
determine the solution, for example the net flux and flux density of an initial, homoge- 
neous vertical magnetic field. Boundary conditions, therefore, remain a hot topic also in 
future. 

Most excitement in carrying out numerical simulations comes from the prospect of 
performing experiments with the object under investigation: experiments in the numeri- 
cal laboratory. Not only that an astrophysical object can be reconstructed and simulated 
in the virtual world of the computational astrophysicist. Once in the computer, he has 
the prospect of carrying out experiments as if the celestial body was brought to the lab- 
oratory. 

The following few examples shall illustrate some aspects of this. 



More detailed physics 

In the solar chromosphere the assumption of LTE (local thermodynamic equilibrium) 
is not valid. Even the assumption of statistical equilibrium in the rate equations is not 
valid because the relaxation time-scale for the ionization of hydrogen approaches and 
surpasses the hydrodynamical time scale in the chromosphere [IFTll . Yet, in order to take 
time dependent hydrogen ionization in a three-dimensional simulation into account, 
simplifications are needed. Leenaarts and Wedemeyer-Bohm Il68ll employ the method 
of fixed radiative rates for a hydrogen model atom with six energy levels in the three- 
dimensional radiation (magneto-)hydrodynamics code CO^BOLD. Thus, additional to 



the hydrodynamic equations, they solve the time-dependent rate equations 
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with Pjj being the sum of coUisional and radiative rate coefficients, Pij = Cjj+Rij. The 
rate coefficients are now local quantities given a fixed radiation field for each transition, 
which is obtained from one-dimensional test calculations. 

Simulations with this approach show that above the height of the classical temperature 
minimum, the non-equilibrium ionization degree is fairly constant over time and space 
at a value set by hot propagating shock waves. This is in sharp contrast to results with 
LTE, where the ionization degree varies by more than 20 orders of magnitude between 
hot gas immediately behind the shock front and cool regions further away. The addition 
of a hydrogen model atom provides realistic values for hydrogen ionization degree and 
electron density, needed for detailed radiative transfer diagnostics. Using this method, 
Wedemeyer-Bohm et al. [|69ll have computed synthetic maps of the continuum intensity 
at wavelengths around 1 mm for a prediction of future high resolution maps that are 
expected to be obtained with the Atacama Large Millimeter Array (ALMA). 



Large box simulations 

Benson et al. [70| have carried out first simulations with a large simulation box so as 
to accommodate a supergranulation cell. They started a simulation that encompasses a 
volume of 48 x 48 x 20 Mm-^ using 500-^ grid cells. With this simulation they hope to 
find out more about the large scale motions in the solar convection zone and to carry out 
helioseismic experiments [TTTI 17211 . Ustyugov (46] carries out a similar calculation in a 
box of (18 Mm)^ for which he uses 192 x 192 x 144 grid points. 

Hansteen [73] has carried out MHD simulations comprising a vast height range from 
the top layers of the convection zone into the transition region and the corona. With these 
simulations they seek to investigate various chromospheric features such as dynamic 
fibrils [74.1 , mottles, and spicules, which are some of the most important, but also most 
poorly understood, phenomena of the Sun's magnetized outer atmosphere. But also 
the transition zone and coronal heating mechanisms are in the focus of these kinds of 
"holistic" simulations. 



Improvements in boundary conditions 

Many conventional magnetohydrodynamic simulations of the small-scale solar mag- 
netic field assume that the horizontal component of the magnetic field vanishes at the 
top and bottom of the computational domain [e.g.|75l|76l|39ll4llj which is a rather stark 
constraint, especially with respect to magnetoacoustic wave propagation and Poynting 
flux. Recently, Stein and Nordlund [42J have introduced an alternative condition with 
the possibility of advecting magnetic field across the bottom boundary (see also Fig. |6] 
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FIGURE 16. Left: Still image of the logarithmic magnetic flux density from a time series for the 
instant t = 1368 s after starting with an initially homogeneous vertical field of 1 mT flux density. A strong 
magnetic flux sheet has formed extending from (x,z) « (2000, —500) to {x,z) ~ (2500,0). Right: A plane 
parallel wave with frequency 100 mHz travels through convecting plasma into the magnetically structured 
photosphere and further into the low )3 (magnetically dominated) chromosphere. The panel shows the 
difference in absolute velocity between the perturbed and the unperturbed solution 212 s after the start of 
the perturbation. The wave becomes strongly refracted in the low )3 region and at the location of the flux 
sheet. From Steiner llT4l . 



right). Thus, upflows into the computational domain carry horizontal magnetic field of 
a prescribed flux density with them, while outflowing plasma carries whatever magnetic 
field it instantly has. With this condition an equilibrium in which equal amounts of mag- 
netic flux are transported in and out of the computational domain is approached after 
some time. It should more faithfully model the plasma flow across the lower boundary 
and it allows for the effect of magnetic pumping 11771 1781 . 



Helioseismic experiment with a magnetically structured atmosphere 

With numerical experiments Steiner et al. [|79l have explored the feasibility of using 
high frequency waves for probing the magnetic fields in the photosphere and the chromo- 
sphere of the Sun. They track an artificially excited, plane-parallel, monochromatic wave 
that propagates through a non- stationary, realistic atmosphere, from the convection-zone 
through the photosphere into the magnetically dominated chromosphere, where it gets 
refracted and reflected. 

When comparing the wave travel time between two fixed geometrical height levels 
in the atmosphere (representing the formation height of two spectral lines) with the to- 
pography of the surface of equal magnetic and thermal energy density (the magnetic 
canopy or /3 = 1 surface) they find good correspondence between the two. These numer- 
ical experiments support expectations of Finsterle et al. [80J that high frequency waves 
bear information on the topography of the 'magnetic canopy' . This simulation exempli- 
fies how a piece of Sun can be made accessible to virtual experimenting by means of 
realistic numerical simulation. 



PART II: THEORETICAL ASPECTS OF MAGNETIC FLUX 

TUBES 



In this second part we discuss some basic properties of small-scale magnetic flux tubes. 
In order to do so, we abstract from the complexity and peculiarities of observed magnetic 
flux concentration or of simulations as the ones shown in the first part of these lecture 
notes. We will consider ideal flux tubes with a sharp boundary and often with a circular 
cross section. Although such ideal flux tubes do not exist in the non-ideal plasma of the 
solar photosphere, the approximation is good enough for this idealization to highlight 
some basic properties that also pertain to more realistic magnetic flux concentrations. 

After definition of an ideal magnetic flux tube we will compute the current sheet that 
is associated with the tube surface and will construct a hydrostatic, vertical flux tube and 
study thermodynamic properties of it. We will derive a criterion for stability with respect 
to the interchange or fluting. 

Supplementary chapters of the lectures that go beyond the scope of these notes on 
flows in a magnetic flux tube and how tubes of superequipartition field-strength can be 
formed and on wave modes supported by magnetic flux tubes can be found, e.g., in the 
book of Stix MB 



CONCEPT AND PROPERTIES OF MAGNETIC FLUX TUBES 



A magnetic flux tube or magnetic flux bundle is defined by the surface generated by 
the set of field lines that intersect a simple closed curve. The magnetic flux, crossing a 
section S of the flux tube is given by 



From the divergence free nature of the magnetic field and using Gauss's law we obtain 




(14) 




Bds+ / Bds-F / B-nda = 0, 
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(15) 




FIGURE 17. Segment of a magnetic flux tube. After Priest ll82l . 




FIGURE 18. Coronal loops over the solar limb on November 6, 1999 taken in the 17.1 nm passband 
with the Transition Region And Coronal Explorer (TRACE). Cortesy, Stanford-Lokheed Institute for 
Space Research and NASA. 



where we use that the flux tube surface is by definition parallel to the magnetic lines of 
force, hence, B ■ n = 0. As a consequence of V • B = 0: 

Fi=F2. (16) 

Hence, like a water hose that keeps mass flux constant, a magnetic flux tube conserves 
the magnetic flux. 

Flux tubes can be considered to be the building blocks of solar magnetism. However, 
they must not necessarily be thought of as independent isolated structures. Coronal loops 
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FIGURE 19. Pore (left) and corresponding Dopplergram (right). A downdraft (red) exists at the periph- 
ery of the pore. See internet version for colors. Courtesy, Kashia Mikurda, Kiepenheuer-Institut. 



for example, which are close to ideal flux tubes, are embedded in an environment that 
may harbour magnetic fields of similar strength. They just appear like "isolated" flux 
tubes because the matter within it happens to be hot enough to emit radiation to become 
visible in certain optical passbands. Figure [TS] shows examples of coronal magnetic flux 
tubes over the solar limb seen face-on and edge-on. 

Pores are another good example of magnetic flux tubes. They often have a circular like 
cross section and seem to extend in close to radial direction from the convection zone 
into the photosphere. Unlike sunspots pores show no penumbra. Fig. [19] (left) shows 
a pore with granulation immediately bordering its periphery, giving the impression of 
an "isolated" magnetic flux tube embedded in a field-free or weak field environment. 
Fig. 19 (right) is the corresponding Dopplergram, showing a downdraft in the periphery 
of the pore very much like the downdraft indicated in Fig. [2j 

Real magnetic flux tubes on the Sun are probably not ideal in the sense of Eq. (16) 
but are more likely "leaky" ([83J), viz., magnetic lines of force do leave and enter 
the tube keeping Eq. (16) only approximately conserved. Also from multidimensional 
simulations we know (see chapter on results from multidimensional simulations) that 
often field lines of a flux concentrations near the surface disperse again in shallow 
depths so that real flux tubes can be expected to have "open ends". Hence, magnetic 
flux concentrations at the solar surface may have a flux tube character over only a short 
distance, typically around the level of optical depth unity. 



What confines a magnetic flux tube? 

In the following we take a closer look to the flux-tube surface, where a discontinuity 
in magnetic field strength and gas pressure exists. 

Consider a magnetic flux tube, embedded in a field-free or weak-field medium. The 
discontinuity in magnetic field strength at the flux-tube boundary can be formally de- 
scribed with the step function 6{x) with 6{x) =0 for x < and d{x) = I for x> O.ln 
the coordinates of the local frame given by s and n, where s is tangential to the magnetic 
field of the flux-tube surface (see Fig. |20|), B is given by B = (0, 0, 5,- - [5, - 5^] (^ ) ) |^ 



Applying Ampere's law: 

VxB = Aij (17) 
and using d'{£,) = 5{^) (Dirac's 5-distribution) we get: 

j = i(0,[5]5(^),0), (18) 

where [B] = Bi — Bg. Integration over an £-range and letting £ ^ 0, leads to the current 
sheet which flows perpendicular to the s-n plane: 

r = -[5]. (19) 



A more refined formulation with an expansion of B, and around the origin of the local frame yields 
identical results. 



flux tube 




FIGURE 20. Flux-tube boundary with local coordinates. § is tangential to a field line on the tube surface, 
n perpendicular to it (main normal). The field strength at the tube surface immediately inside and outside 
of it are Bj and Be, respectively. 



Notice that this is not a regular current density since we have carried out a spatial 
integration in the coordinate ^ in the direction of n across the flux-tube boundary so 
that j* has dimension Am^^. 

Thus, we have obtained as a first result that a current sheet occurs at the flux-tube 



surface with the current flowing perpendicular to the magnetic lines of force. Fig. 15 
nicely demonstrates how current sheets may occur in the actual solar atmosphere. The 
characteristic double layers in the panel corresponding to 90 km above the mean surface 
of optical depth unity originate from the current sheet that flows on either side (surface 
with a sharp gradient in B) of the magnetic flux sheets that are located in the intergranular 
lanes. 

Next, we consider the force balance starting from the momentum equation 

Dv 

P^ = -Vp + (jxB)+pg, (20) 

where the middle term on the right hand side is the Lorentz force and g the gravitational 
acceleration at the solar surface (optical depth Tsoonm = !)• Again using Ampere's law. 



Eq. ( 17 1, we obtain for the static case 



Vp = -((VxB)xB)+pg (21) 
and further using the vector identity (V x B) x B = (B ■ V)B - ^V(B ■ B), 

V/' = -V(^)+pg+i(B.V)B. (22) 

Now we decompose the last term into a component parallel (s) and perpendicular (main 
normal n) to a surface field line. From B = 5s follows: B ■ V = 5 • s ■ V = B{d / ds) and 
further 

(B.V)B = (5^)B = 5^ (5s) = A + = li^)s + B'^ , (23) 

OS OS OS OS OS 2 Rr 



where Rc is the curvature radius of the field line. Thus, we obtain: 



Vp=-V(-) + pg+^(-)s+--, (24) 



where we note that the term B^/lji acts like a gas pressure, called the magnetic pressure 
[see im chap. 2.7]. Multiphcation of this equation with n yields the force balance 
perpendicular to the tube surface: 

dp d B^^ . 52 1 

and integration over a small interval [— £, e] across the tube surface gives: 

5?/2M 2 -2 




■ ^ 



-^Bl-B]) 

Letting e ^ 0, we finally obtain: 

P. + ^=P, + ^. (26) 

Thus, we have obtain as a second result that a discontinuity in pressure occurs at the flux- 
tube surface so that the Lorentz force is balanced by a corresponding pressure force. It is 
important to realize that Eq. ([26]) was obtained from basic principles without using the 



"thin flux-tube approximation". Eq. (26) expresses that the total pressure consisting of 
gas pressure plus magnetic pressure is constant across the flux tube surface. This results 
holds generally for tangential magnetic discontinuities in a plasma. Note however, that 
indices / and e refer to locations in the immediate vicinity of the discontinuity. 



A microscopic picture of the current sheet 

It may seem surprising that an electrical current can be confined to a thin sheet in a 
bulk of plasma with high conductivity everywhere in space. This subsection is intended 
to sketch the microscopic processes that lead to this current sheet. Consider first the 
situation in which we have a uniformly directed magnetic field with a spatial gradient 



in field strength as plotted in Fig. 21 (left). The field is directed perpendicular to the 
drawing plane pointing out of it and increasing in strength from bottom to top. A 
negatively charged particle that gyrates in this field will unabatedly drift to the right 
because the Lorentz force, = q{\ x B), is not constant along the path of the particle 
but increases as it moves into higher field strength leading to a reduction in curvature 
radius. Likewise, a positively charged particle will drift to the left. Hence, electrons and 
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FIGURE 21. Left: Particle motion in a magnetic field with a transverse gradient showing grad B drift 
motion (adapted from Krall and Trivelpiece |85 1). Right: Particle orbits in a plasma with a sharp density 
gradient resulting in a net drift current (right). 



ions in a magnetic field with a transverse gradient in flux density show a gradient B drift 
motion, leading to a net electrical current. 

Next we consider in Fig. |2T] (right) the situation of a uniform magnetic field pointing 
out of the drawing plane. Instead of the gradient in field strength we have a contact 
discontinuity in density with a lower density in the lower half of the figure and a 
higher density in the upper half. In each half, positively charged particle will gyrate 
in a circle in the sense like indicated in Fig. |2T| (right), however, there are much more 
gyrating particles in the upper half than in the lower half because of the difference 
in particle density. Considering now the hatched control volume, which straddles the 
contact discontinuity and which is much smaller that the gyration radius, we can readily 
see that much more charged particles run to the left than to the right, leading to a net 
positively charged current to the left and, correspondingly, a net negatively charged 
current to the right. Due to the imbalance of gyrating particles, an electric current results 
without a net transport of charges. No net current results for a control volume in either 
the upper or the lower half if it does not include the discontinuity. 

Finally, Fig. |22] sketches the situation at the boundary of a flux tube that is thought to 
be embedded in a field free environment. Charged particles in the field-free plasma move 
on straight free paths but start to gyrate as soon as they enter the magnetic field of the flux 
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FIGURE 22. Particle orbits in a plasma with a sharp gradient in density and magnetic field strength 
as it occurs on the surface of a magnetic flux tube. It results a net drift current confined to the flux-tube 
surface. 



tube, where negatively charged particles will drift to the right for the situation depicted 
in Fig. 22 However, this drift current will not be cancelled by gyrating particles within 
the flux tube because of the reduced particle number-density there. Hence, we obtain a 
net drift current with a net transport of charges. 



THE EQUATIONS FOR A HYDROSTATIC FLUX TUBE 



In this section we derive the equations that govern the magnetic structure of a vertical, 
axisymmetric, untwisted magnetic flux tube embedded in a gravitationally stratified, 
field-free atmosphere. Although such a perfectly symmetric object is an idealization 
that does probably not occur in reality, the basic properties of the solutions to these 
equations can be recovered in the more complex situations that we know from numerical 
simulations in a three-dimensional environment. The system of the magnetohydrostatic 
equations to be solved is 
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-V;? + pg+jxB 
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(27) 
(28) 
(29) 



We decompose Eq. (27 ) in components parallel and perpendicular to the magnetic field. 
Thus, multiplying Eq. (271) with B gives 



B-(Vp-pg)=0 



(30) 



and taking the cross product of Eq. ( 27 ) with B gives 

. 1 



52 



Bx(Vp-pg) 



(31) 



using that B x (j x B) = B^j — (B ■ j) ■ B = B^j for an untwisted axisymmetric field. 



Next we reform Ampere's law, Eq. ( [28] ). Using cylindrical coordinates and Eq. ( |29| ) 
we first have 



B = Vx A: 



which reduces to 
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(32) 



(33) 



for an axisymmetric flux tube without twist. With the stream function ^ 
write 
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r dz 
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and further, taking the curl of B 
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FIGURE 23. Magnetic field line ^ with local coordinate frame, s is tangential to the field line, n is the 
main normal. 



SO that Amperes law for the (^-component (the r- and z-components being zero), be- 
comes (Grad-Shafranov): 



1 ^2vj/ 

r dr dz^ 



(36) 



This is a second order elliptic, inhomogeneous, partial differential equation. It is non- 
linear because depends on ^ again as we will see further down in the text. The 
magnetic field components can be recovered from the stream function ^: 
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r dz 
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r dr 



(37) 



Note that contours of ^ are field lines of the system. 



Now we consider Eq. (30) and Fig. 23 If gravity acts parallel to the z-axis in the 
negative direction and s measures the distance along a magnetic field line inclined at an 



angle to the vertical direction, we get from Eq. (30) 



dp 

Ts 



+ pgcos6 = 0. 



Substituting the ideal gas equation for the density p we obtain 



dp 

P 



mg 



dz. 



(38) 



(39) 



H = kgT /{mg) is a length scale, called the pressure scale height, because the gas 
pressure drops by a factor e over this distance in the outward direction. Note that H 
is proportional to the temperature, m is the mean molecular weight. Integration from a 
reference height z = to an arbitrary height z along a magnetic line of force labeled with 
^ gives 

dz' 



7? = Po(^)exp 



H{T{^>,z')) 



(40) 



This equation expresses that hydrostatic equilibrium holds along magnetic field lines. 
Picking a field line, we obtain with this equation the gas pressure, p{z = 0), exerted by 
the mass column along this field line. 

Having made use of the longitudinal component of Eq. ( |27| ) we now turn to its 
perpendicular component, Eq. (31 1. Using cylindrical coordinates B = {Br , , By) and 
Vp = {dp/dr , , dp/dz), Eq. (31 1 becomes 



1 

52 



B 



dp 



(41) 



From Eq. (40 1 we obtain 
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and 
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dr 



(42) 



so that the above equation for reduces to 



dp^ 
(9^ 



(43) 



At the surface of the tube the magnetic field will, in general, be discontinuous, resulting 
in a sheet current there. Using 



1 



Pe-Pi = :^iBt-Bi) 



1 



[B,-Be){B,+Be) 



and Eq. ( 19 1 we finally get 



^{Pe-Pi) 



(44) 



(45) 



B, + Be 

The procedure for computing the magnetic structure of a vertically oriented, axisym- 
metric flux tube without twist is the following: 

(i) Specify an initial magnetic configuration, e.g., with the help of the thin flux-tube 
approximation. This initial flux tube shall be embedded in a given plane-parallel 
atmosphere that is thought to remain unchanged. 

(ii) Calculate the gas pressure everywhere in space using hydrostatics along magnetic 
lines of force, Eq. ( [40] ). Since we do not solve an energy equation here, the pressure 
scale height as a function of field line and height, Hi^.z), must also be specified. 
An acceptable guess would be that H is independent of ^ and identical to the scale 
height in the surrounding field-free atmosphere, which is approximately fulfilled if 
the temperature is constant in planes parallel to the solar surface. The pressure and 
pressure scale height in the surrounding field-free atmosphere is known from the 
initial, embedding atmosphere. 



(iii) Evaluate volume and current sheet using Eqs. ( [43[ ) and ( [45| ). Eq. ( [43[ ) applies to 
regions of continuous gas pressure, Eq. ( |45| ) to locations of discontinuous gas 
pressure and field strength like the surface of the magnetic flux tube. 
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FIGURE 24. Boundary conditions for solving the elliptic partial differential equation, Eq. ( 36 1 



(iv) Solve the Grad-Shafranov equation, Eq. ( |36| ), with appropriate boundary condi- 
tions. This leads to a new field configuration, ^(r,z), that must not necessarily be 
identical to the initial magnetic configuration. Hence, one must go back to (ii) and 
iterate (ii) to (iv) until convergence has been achieved. 



Boundary conditions need to be specified for solving Eq. (36). A plausible choice 
is shown in Fig. 24 The field line with ^ = that marks the left boundary of the 
computational domain coincides with the symmetry axis of the flux tube. At the bottom 
boundary the vertical component of the magnetic field as a function of radial distance 
is prescribed, which fixes ^(r,z = 0). The boundary on the right hand side is thought 
to separate the flux tube from fields of neighbouring flux tubes. We can think of it as a 
separating field line of constant ^. At the top we assume that the magnetic field merges 
with the magnetic field of similar neighbouring flux tubes leading into a homogeneous, 
vertical magnetic field with vanishing horizontal component. Formally, this condition 
is obtained by setting d^'/dz = —rBr = 0. The width W of the computational domain 
determines the filling factor /, viz., the fraction of the area at z = that is occupied with 
magnetic field: 

f = W^/Rl. (46) 



Eq. (36) together with these boundary conditions constitutes a free boundary problem 
because the boundary of the magnetic flux tube is not known from the beginning but is 
part of the solution. 

Instead of using a rectangular computational domain, Pizzo ll86l uses a "body-fitted" 
nonorthogonal coordinate system to map the physical domain of the flux tube into a 
unit square computational domain. A multigrid elliptic solver is used at each iteration 
stage for solving the Grad-Shafranov equation. Fiedler and Cally [|87l use a similar mesh 
in which contours of constant ^ (field lines) constitute one coordinate, the normalized 
arc length along field lines the second one. Jahn [88] and Jahn and Schmidt [i89il use a 
similar method for the construction of sunspot models. 



In case of horizontal temperature equilibrium, 7^ (z) = Te{z) = T{z), we have (neglect- 
ing any horizontal variation in ionization degree) Hi{z) = He{z) = H{z). Then, 
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Pe{z) =POeexp 
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(47) 



from which follows that piiz) < Pe{z)'^z assuming that pi does not depend on radius 
(thin tube approximation) because poe ~ PQi = Bq/2ii > 0. Since Ti{z) = Te(z) it follows 
for the densities that p,(z) < Pe(z) Vz. Since above conditions can be expected to hold 
very well in the photosphere we can say that photospheric flux tubes are rarefied, one 
also says "partially evacuated". 

As a most simple case we explore in the following a flux tube with a constant axial 
field strength and gas pressure at the base level, zo, where the radius is Rq. Then 
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from Eq. (37) and 



PiO = PeO — 



bIo +bIq{Ro) 
111 



(48) 



(49) 



from the requirement of pressure balance, Eq. (26 1. We note, that the final B,-q{Rq) is not 
known from the very beginning so that the boundary condition for the pressure needs to 
be adjusted in the course of the iteration. 

If we further assume that the temperature at a given hight level is constant, then the 
gas pressure is constant too, so that the volume current is 



dp 



= 



(50) 



and we have a potential field inside the flux tube. However, there remains a current sheet 
at the tube surface: 



Fig. 25 shows the corresponding solution for a flux tube with a field strength of 0.15 T 
and a radius of 100 km at the base of the photosphere (Tsoonm = !)• The figure shows 
representative lines of force of the magnetic field. We recall that the left hand boundary 
coincides with the symmetry axis of the cylindrical tube. The ring areas between two 
field lines contain equal amounts of magnetic flux. 

Superimposed on the field-line plot are plots of the radial variation of Br (dashed 
curve) and of (solid curve) at different heights, both normalized to the value of B^ at 
the axis, indicated in gauss (10^ T). The flux tube merges with the field of neighbouring 
flux tubes at a height of ^ 500 km. The filling factor is / = {Rq/W)^ = 0.1. 

Due to the exponentially decreasing gas pressure in the gravitationally stratified at- 
mosphere, the magnetic flux tube must expand with height because of pressure balance, 
Eq. ( [26] ), and flux conservation resulting from VB = 0. Even though, the axial field 




FIGURE 25. Field lines and cross sections for a cylindrically symmetric, vertical flux tube with uniform 
pressure and axial field strength at the base, carrying a current sheet at the surface. Superimposed on the 
figure are plots of the radial variation of B, and at different heights, normalized to the value of B^ at the 
axis, each indicated in gauss (10^ T). Note that the vertical scale is compressed relative to the horizontal 
scale by a factor 1.6. From Steiner et al. |,90J . 



Strength is constant with radius at the reference height z = 0, this is no more the case 
higher up in the photosphere, where By{r) decreases towards the flux-tube boundary. 
This variation of the axial field strength with flux-tube radius may become much stronger 
than in the case of Fig. |25] depending on the expansion rate of the flux tube, which in 
turn depends on the combination of internal to external atmosphere. If Hi = He as was 
assumed for the flux tube of Fig. [25] it can be shown with the help of the thin flux-tube 
approximation that it takes about four pressure scale heights for the flux-tube radius to 
expand by a factor of e and about two pressure scale heights for the axial field strength 
to decrease by a factor of e. 

Introducing the scalar function G = one can treat the case of a twisted axisym- 
metric flux tube. The Grad-Shafranov equation then becomes 

fdG ^2vj; i^vj/ ^2vj/ 

\ dz ' dr-^ r dr dz^ ' dr 




If G = G(^) (torque-free condition) we need only solve the ^-component, for which: 



In the absence of an external magnetic field, the magnitude of the current sheet remains 
the same as before: |j*| = 2{pe — pi)/Bi, directed perpendicular to the field lines at the 
surface and, hence, no longer purely azimuthal. The ^-component of the current sheet 
is: 

The solution procedure remains the same as before. Without even carrying out a stability 
analysis it can be shown that there exists an upper limit for twisting the flux tube beyond 
which, there is no magnetohydrostatics possible because the magnetic tension at the 
tube surface can no longer be counterbalanced by gas and magnetic pressure in the tube 
interior. 



The magnetic canopy 



The expansion rate of hydrostatic flux tubes critically depends on the combination 
of external to internal atmosphere. This is demonstrated by the following ideal model 
in which we assume constant but different pressure scale heights, Hg and //, in the 
atmospheres external and internal to the magnetic flux tube, respectively. If Hi > He, 
which is equivalent to a hot flux-tube atmosphere, 7^ > Tg, a critical hight, Zcrit^ ensues, 
where the flux tube must strongly expand in the horizontal direction. 
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FIGURE 26. Gas pressure as a function of height within (i) and external (e) to the flux tube. At the 
critical height, Zcrit^ the magnetic field is forced to expand in the horizontal direction. 



This situation is illustrated in Fig. 26 which shows the logarithmic gas pressure 
normalized to the external gas pressure at the reference height, zo^ as a function of 
height in the atmosphere. Initially, the internal gas pressure, pi, is smaller than the gas 
pressure in the external atmosphere, pe, for to accommodate the magnetic flux tube with 
magnetic pressure pmag = Pe — Pi- But at the critical height, Zcrit, the internal gas pressure 
approaches the external gas pressure so that the magnetic field (magnetic pressure) at the 
tube surface must vanish, which causes the sudden horizontal expansion. 

Generally, the pressure stratification in the internal and external atmosphere are given 

by 

Pi,e{z) = POi,ee^P\ - J JfJ^ I • (55) 
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FIGURE 27. Flux tube with a standard quiet sun atmosphere embedded in a radiative equilibrium 
atmosphere without a chromospheric temperature rise. The field-strength at the base (z = 0) is 0.15 T. 
The field lines spread into a horizontally extending canopy field at a height of 900 km. 



With Peizcnt) = Piizcnt) and assuming constant pressure scale heights, we obtain 
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In = — — In — - — , (56) 



Hi — He poi Hj — H, 



where we use the plasma beta, /3o 
the reference height zq- 



POi/ Pmag, the ration of gas to magnetic pressure at 



As an example. Fig. 27 shows magnetic lines of force of a rotationally symmetric, 
vertical flux tube that harbours a standard quiet sun atmosphere (C' of Maltby et al. [|9T| , 
Table 11). The external atmosphere is an atmosphere in radiative equilibrium showing 
no chromospheric temperature rise [92]. Therefore, above the temperature minimum of 
the internal atmosphere, at a height of approximately z = 500 km, it becomes distinctly 
cooler than the flux-tube atmosphere with the consequence that p, rapidly approaches 
Pe so that the field lines spread into the horizontal direction at a height of about 900 km 
forming a so called canopy field. The strength of the canopy field is weak - in fact, the 
magnetic pressure must vanish at the boundary to the canopy as we have seen. Therefore, 
the plasma /3 is very inhomogeneous within a tube like that of Fig. 27 and the surface 
of /3 = 1 spreads much less dramatic with height than the peripheral field lines might 
suggest. 

Canopy fields also occur in the superpenumbra of sunspots or in the surroundings of 
network magnetic fields. The latter can be directly observed in limb magnetograms of 



chromospheric spectral lines. This is illustrated in the sketch of Fig. 28, which shows an 
expanding magnetic flux tube that should be thought to represent a network magnetic 
flux concentration. The solar surface (T500 = 1) is inclined with respect to the horizontal 
to indicate observation near the solar limb. The hatched region is filled with magnetic 
field of the flux concentration, while indicates the field-free region. 




FIGURE 28. Sketch of a chromospheric magnetogram near the solar limb (a) and the corresponding 
magnetic configuration (b). The base of the magnetic canopy is located at a height Zc above the solar 
surface. From Steiner [93 J . 



Lines of sight passing through the central network field concentration receive a mag- 
netic field component toward the observer and so do lines of sight passing through the 
canopy field to the left (disk center) side of the network (dashed line of sight). The 
canopy field to the right (limbward) side of the network, however, gives rise to a line 
of sight component of opposite direction (solid line of sight). Thus, the observed mag- 
netic polarity inverses across a line that coincides with the limbward edge of a unipolar 
magnetic network field as is sketched in the corresponding magnetogram to the left. 

Indeed it was observed that chromospheric magnetograms of a network region that is 
unipolar in magnetograms of photospheric lines (that form below the canopy) show near 
the limb a fringe pattern in polarity [94|. This fringe pattern can be easily understood 
when we imagine the sketch of Fig. |28] (right) to be continued to a sequence of neigh- 
bouring unipolar flux concentrations, each giving rise to a chromospheric magnetogram 
like the one of Fig. [28] (left). From magnetograms taken in the line of Mg Ib2 (517.3 nm) 
Jones and Giovanelli [|95l derived canopy heights of 600-1000 km in quiet-Sun regions 
and as low as 200 km in active regions. 



MAGNETIC FLUX TUBE IN RADIATIVE EQUILIBRIUM 

There are two basic modes of energy transport in the solar photosphere and convection 
zone: radiative and convective. When in a stationary state all of the energy is transported 
by radiation, we have radiative equilibrium, conversely, in the case of pure convective 
transport we have convective equilibrium. In a stationary transport process, the frequency 
distribution of the radiation, or the portioning of energy between the radiative and the 
convective mode of transfer, may be altered; but the energy flux as a whole is rigorously 
conserved. Formally, this is expressed by 



V-Ftot = 0, where Ftot = Frad + F, 



(57) 



assuming that no energy transport other than radiative and convective be relevant, as is 
usually assumed in mixing length theory. 

For the two limiting cases mentioned above we have 



radiative equilibrium: V ■ F^ad = Fcon = , (58) 
convective equilibrium: V ■ Fcon = Frad = . (59) 

In the solar photosphere radiative energy transfer by large prevails so that V ■ Frad = 
is a good approximation. With /(r, n, v) being the radiative intensity with frequency v 
propagating in direction n at location r in three-dimensional space, the total radiative 
flux is given by 

' '^/(r,n,v)ndvd(0. (60) 



rad ~ 

Ak 

/(r, ft, v) has dimension Wm^^hz^^sr^^ and correspondingly Fj-ad has dimension 
Wm^. The radiation field follows as a solution of the radiative transfer equation in 
three-dimensional space [|96l 

(ft ■ V)/(r, ft, V) = T] (r, V) - K(r, v)/(r, ft, v) . (61) 

The emissivity r/(r, v) is given by T](r, v) = K"(r, v)5(r, v) with S being the source 



function. If 5(r) and K'(r) are known, Eq. (61) is an ordinary first order differential 
equation whose solution is called the formal solution of the radiative transfer equation. 
From these equations we obtain: 

oo 

V-Frad = yy (ft ■V)/(r,n,v)dvdw (62) 

o 

(K(r, v)5(r, v) - K(r, v)/(r, ft, v) = . (63) 
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/(r,v) = ^y'/(r,ft,v)d£0 (64) 

An 

being the mean intensity, we finally obtain the constraint equation for radiative equilib- 
rium: 

oo oo 

j K:(r,v)5(r,v)dv = j K:(r, v)7(r, v)dv . (65) 



In general we do not know the temperature distribution r(r) that satisfies radiative 
equilibrium. If we start with a guess r(*^)(r) for which we have calculated the correct 
source function 5(r, v), we will find that the constraint equation for radiative equilibrium 
is not satisfied. It is therefore necessary to iteratively adjust r(r) until the requirement 
of radiative balance is satisfied. 




FIGURE 29. Two-dimensionally structured atmosphere in radiative equilibrium. Left: Radiative bound- 
ary conditions (for determining G(r) of Eq. 66 1 and distribution of opacity. Right: Isotherms of the 
radiative equihbrium solution. From Steiner 11971 . 



Assume a given temperature distribution r(r) of the magnetohydro static configura- 
tion for which we also know the pressure p(r) and density p(r). This allows us to derive 
LTE- values for the opacities K"(r, v). We then are able to compute the radiation field 
everywhere in space r and for all frequencies v by evaluation of the formal solution of 
the radiative transfer equation 

7v(r) =Av(r,r')5v(rO+Gv(r) , (66) 

where Ay is the integral operator which adds the intensities at r caused by emission at all 
the points r' in the considered computational domain, and where Gy is the transmitted 
mean intensity due to the incident radiation field into this domain. We also use strict 
local thermodynamic equilibrium (LTE) (5v = By), although a scattering component 
could be included in the solution method that follows. Note that Ay together with Gy 



is just a short-hand for the integration, Eq. 64, and the formal solution of the radiative 



transfer equation, Eq. (61 ). In practice we can think of Ay as a matrix operator acting on 
a vector of source values given at discrete spatial locations. 

Defining the integral operator ^ so that ^0 := /J°K'y0dv , where ^ is a scalar 



function, we can write the constraint equation for radiative equilibrium, Eq. (65 1, as 



jXTBy = JTAyBy + J^Gy . (67) 

Given an initial temperature distribution, T^*^) (r) , this equation is in general not satisfied. 
But we can compute a correction Ar(r) such that 

jr5y(r('^)+Ar) = jrAy5y(r(°)) + jrGy (68) 

is satisfied. Expanding By{T + AT) ^ By{T) + {dBy/dT)\TAT, we can compute the 
temperature correction: 

J^iAy-l)By{T)+J^Gy 

""^^ Jf{dBy/dT)\r • ^^^^ 
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FIGURE 30. Radiative equilibrium solution for a cylindrically symmetric, vertical magnetic flux tube 
with a field strength of sa 0. 15 T at z = 0. The horizontally running contours are isotherms labeled with the 
temperature in K. The curves labeled in the right margin of the figure are contours of constant optical depth 
for vertical lines of sight. The temperature in the convectively unstable region of the indicated pentagonal 
torus was prescribed. Adapted from Steiner and Stenflo |98j|. 



For most practical purposes this iterative scheme is very slow. The next better scheme 
to use is a Jacobi-like iteration, in radiative transfer known as accelerated A-iteration : 

^ jr(Av-i)^v(r) + jrGv 

jf{l-X*){dBy/dT)\T ' ^ ^ 

where Xy is the diagonal element of the matrix representing Av(r). Note that AT — 
Ar(r) and that the above equation must be evaluated for each location r. Since the 
formal solution of the transfer equation must be computed in each iteration at each 
location r, an efficient solver is needed, which is given by the so called method of 
short characteristics [|99l II 0011 and by exploiting the symmetry properties in the case 
of a rotationally symmetric flux tube HlOlll . 



Fig. 29 is a striking example of a two-dimensionally structured atmosphere in radia- 



tive equilibrium. The panel to the left shows the boundary condition (for determining 



G(r) of Eq. 66 ) and the distribution of opacity. The rectangle in the upper left comer 
may represent a magnetic flux sheet with a rarified atmosphere and correspondingly 
reduced opacity. The incident radiation on the bottom and on the right side of the com- 
putational domain is computed from the undisturbed, plane-parallel atmosphere with 
K= Ke= const. The left side is an axis of symmetry of the configuration. The top bound- 
ary has no incident radiation. The panel to the right shows the isotherms of the radiative 
equilibrium solution obtained with the above given solution procedure. While in deep 



P4 



FIGURE 31. Irradiation on point Pi located on the flux-tube axis above the level of Tc = 1 and on point 
P2 located at the same geometrical height but laterally separated from the flux tube in the undisturbed 
atmosphere. Adapted from Steiner II 171 . 



layers the rarified region is cooler than it is in the unperturbed surrounding atmosphere, 
the top layers are hotter. 

This bahaviour can be recovered when computing the radiative equilibrium solution 
of a magnetic flux tube, like the one shown in Fig. 30 It shows isotherms together with 
contours of constant optical depth for vertical lines of sight, labeled in the right margin 
of the figure. The vertically running curves represent lines of force of the magnetic field. 
The frequency integration in the constraint equation for radiative equilibrium, Eq. ([65]) 
was done using opacity distribution functions of Kurucz [I102II , thus, taking line blan- 
keting effects into account [|98ll . The height scale on the left hand side originates at 
the optical depth Tsoonm = 1 in the surrounding, unperturbed atmosphere, from which 
follows that the computational domain extends into the top surface layers of the con- 
vection zone, where the assumption of radiative equilibrium breaks down. Therefore, 
the temperature structure within the indicated pentagonal torus below z = and outside 
the magnetic field was prescribed and kept constant when computing the radiative equi- 
librium solution. Resort to this stratagem obviously reveals the limits of the radiative 
equilibrium approach. Yet, this solution highlights a typical property of flux-tube atmo- 
spheres that can be recovered in fully self-consistent simulations: the temperature along 
the axis of symmetry increases in the deep photospheric layers less strong with depth 
than in the surrounding, unperturbed atmosphere, hence, the flux-tube atmosphere has 
a "flat temperature gradient". This is a consequence of the "radiative channeling effect" 
by which the photospheric layers of the flux tube get heated through influx of radia- 
tion from the optically denser surroundings. This change in temperature gradient has far 
reaching consequences for the formation of spectral lines within the magnetic flux tube. 
Generally, spectral lines tend to become weaker (line weakening) with the consequence 
that filtergrams of spectral regions with a high density of spectral lines show magnetic 
elements of particularly high contrast (e.g., G-band filtergrams). 

Why is the temperature in the photospheric layers of the magnetic flux tube by about 
200 K higher than in the surrounding atmosphere? Consider two points P\ and P2 located 
at equal geometrical height, P\ on the flux tube axis, P2 far away from the flux tube in 
the unperturbed external atmosphere, as is illustrated in Fig. |3T] Pi receives a higher 
intensity of radiation from the flux tube's hot walls (radiative channeling effect) as 
compared to Pi- Both points receive equal amounts of low radiative intensity from the 
top. Therefore, the mean intensity is higher in Pi than in P2, J\ > J2. Since in radiative 



equilibrium J = S = J B^dv T^, T\ > T2. The opposite behaviour occurs for points 
P3 and ^4 as long as the assumption of radiative equilibrium holds. However, the true 
reason for the reduced temperature in the optically thick layers of the flux tube is the 
inhibition of convective energy transport by the magnetic field. 



INTERCHANGE INSTABILITY OF MAGNETIC FLUX TUBES 



Consider a section perpendicular to the axis of a straight flux-tube as is sketched in 



Fig. 32 A perturbation of the flux-tube boundary in such a way that the volumes V\ and 
V2 are equal, does not change the total energy of the configuration as equal amounts of 
magnetic and thermal energy are exchanged. To zeroth order, the perturbation also does 
not change the energy state of a flux tube embedded in a stratified medium. 
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FIGURE 32. Left: Cross section of an unperturbed cylindrical magnetic flux tube. Right: Same view 
after perturbation showing the initial perturbation of the flutes. Adapted from Krall and Trivelpiece ||85]| . 

The instability that may evolve from this perturbation is called interchange instability, 
as the magnetic field and gas of volume 1 is interchanged with the magnetic field and gas 
of volume 2. It is also called instability because of the shape of the perturbed surface 
(like the vertical parallel grooves on a classical architectural column, called flute). 

Consider a small flux-tube section with a small perturbation ^ (Fig. 33 1, where the 



grey shaded area, V,-, be the magnetic flux tube while Ve is field-free, n is the surface 
normal pointing out of the field-free plasma. In order that there is a net restoring force 
on the displaced surface we must, for the indicated displacement t, , have 



^Pitot > ^Pe 



(71) 



since for the equilibrium configuration po/tot = POe- From this follows the condition for 
stability: 

|^ii|ii-Vp,< |^ii|ii-V(p, + — ). (72) 

This criterion also follows from a more general energy principle due to Bernstein et al. 
[I103II . It is both, necessary and sufficient for stability. 
Within the flux tube we have 



Vp, + V(g)=p.g+i(B.V)B 



(73) 



FIGURE 33. Distorted surface S displaced from the equilibrium position (dashed line). The unit vector 
n is normal to the unperturbed surface and the section is normal to the magnetic field lines close to the 
tube surface. Adapted from Meyer et al. IITSl . 



as was derived in connection with Eq. (22). In the external atmosphere we have hydro- 
static equilibrium 

^Pe = Peg . (74) 



Using these equations in the previously derived stability criterion, Eq. f72| ), we obtain: 

1 



n 



(B-V)B-(p,-p,Og 



From Eqs. ( [73] ) and f74] ) we have 
1 



>0 



(75) 



- [B ■ V)B - (p, - p,)g = V(p, + ^ - P.J , 



(76) 



=0 on s 



which means that the bracketed vector in Eq. (75 ) is parallel to n on the surface S of the 
flux tube. 

Therefore, gravity can be eliminated taking the horizontal component of this vector. 
If h is a horizontal vector pointing out of the flux tube into the field-free plasma we get 



h-[(B-V)B] <0. 



(77) 



This means that along any field line in S the magnitude of the component of B in 
any fixed outward horizontal direction must decrease as a function of height in the 
atmosphere for to have stability. For an untwisted axisymmetric flux tube this criterion 
can be expressed in cylinder coordinates as: 



dBr 



<0. 



(78) 



Fromn-(B-V)B: 



(see Eq. 



23 1, where Rc is the curvature radius of the surface in 



axial direction, results another usefuTTorm of the stability criterion for an axisymmetric 
untwisted flux tube: 

--^ + iPe-Pi)s^^nx>0. (79) 



X is the inclination of S with respect to the vertical direction. 

Using a realistic model atmosphere, Meyer et al. fT5\ came to the conclusion that 
only flux tubes with magnetic flux ^> > 10^^ Wb are stable against the flute instability. 
Sunspots and pores have flux in excess of 10^1 Wb. With typical values of = 100 km 
and 5 = 0. 1 T resulting in ^» ^ 3 ■ 10^ Wb, small-scale magnetic flux concentrations in 
plage and network regions are liable to the flute instability! 

Indeed, from time sequences of G-band-bright-point images (e.g., EIIOl) one gets 
the impression that magnetic elements are subject to fluting (see also the discussion in 
the first chapter of part I). On the other hand their typical life time of 6-8 minutes is 
still in excess of the crossing time for Alfven or sound waves of about 100 s. Also there 
are bright points that undergo continual fragmentation and merging in a relatively stable 
location, persisting over several hours. 
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